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We consider the problem of keeping an arbitrary state ps out of thermal equilibrium. We find that counter¬ 
acting thermalisation using only a resource system which is in a stationary state at the initial time and a system- 
resource interaction that preserves the global energy is possible if and only if the target state ps is block-diagonal 
in the eigenbasis of the system’s Hamiltonian Hs- As a consequence, we compute the extra work the operator 
must provide by tuning the resource-system interaction to overcome this constraint. This quantity, which is 
interpreted as the work needed to preserve the coherences in the state, can be expressed in terms of the target 
state ps and the thermal equilibrium state pp, and it is proportional to the symmetrized relative entropy between 
Ps and pp. 
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-Introduction: The idea of interpreting thermodynamic con¬ 
cepts in the context of quantum mechanics has intrigued the 
physics community for decades. Some of the topics that have 
been proposed in this area of research in the last years include 
quantum applications of Landauer princii jieilHi, the study of 
out-of-equilibrium quantum systems 10,0 the possibility 
of creating quantum thermodynamic engines |0]. Thermody¬ 
namics and quantum mechanics are basically involved in all 
branches of physics, and understanding the interconnections 
between these two theories is of crucial importance for the 
resolution of important problems, ranging from the miniatur¬ 
isation of computer processors 11 to the study of quantum 
behaviours in out-of-equilibrium biological systems (see @] 
for a recent review on various topics in quantum biology). 

On the other hand, in the context of quantum thermodynam¬ 
ics (QT), it is important to determine which states of a given 
system have genuine quantum features. In this regard, the con¬ 
cept of asymmetry with respect to time translation can be used 
to distinguish classical states from quantum statesH. Loosely 
speaking, while it might be possible to interprete the ther¬ 
modynamic properties of states that are symmetric with re¬ 
spect to time translation (i.e. states that are block-diagonal in 
the eigenbasis of the system’s Hamiltonian) using the frame¬ 
work of classical statistical mechanics, the thermodynamics of 
states with non-zero off-diagonal elements (i.e. states with co¬ 
herences) can only be understood in terms of quantum theory. 
Notwithstanding the importance of recent results in the held 
CMl, most of these results focus on the study of block- 
diagonal states, and only very recently the role played by off- 
diag onal elements has become subject of deeper investigations 

lo®. 

Following this line of research, here we consider the prob¬ 
lem of keeping a system 5 in a arbitrary state ps while it is 
in contact with a thermal bath. We assume that the thermal¬ 
isation process is modeled as a thermalising machine, i.e. a 
step process in which each step is characterised by a partial 
swa pping be tween the state of the system and the state of the 
bathi2(X 1^ . In the limit of inhnitesimal steps, this frame¬ 
work reproduces the continuous dynamics of the svstemll2^ 
1^ . We show that the amount of coherence in the state ps 
with respect to the eigenbasis of the system’s Hamiltonian Hg 
crucially determines the properties of the operations needed to 


keep the state out of equilibrium. Specihcally, we show that, 
if we restrict ourself to a certain class of external operations, 
it is possible to keep ps out of equilibrium if and only if ps is 
block-diagonal. We also show that operations keeping states 
with coherences out of equilibrium (thus not belonging to the 
class of operations mentioned above) require an extra amount 
of work in order to be implemented, which can be then inter¬ 
preted as the work needed to keep coherences. 

-Thermalisation processes: We are now going to introduce 
the concept of thermalising machines, which was hrst pro¬ 
posed in 12(1 l2d]l . Here, we will make use of thermalising 
machines in order to give a precise definition of a thermalisa¬ 
tion process. This approach is particularly powerful due to the 
generality of the model, which is applicable to any quantum 
system. 

Let us consider a quantum system S described by a Hilbert 
space T-Ls and with Hamiltonian Hg ■ The system is in contact 
with a thermal bath B at inverse temperature /3. Following 
iSlUl], we describe the Hilbert space of the thermal bath B 
as a collection of N independent copies of the system S. We 
now assume that the initial global state of the system and the 
bath is given by a product state pgp = Ps^pB, with ps an arbi¬ 
trary density matrix for the system and ps = ^k=i 

— (e“^^®)/(Tr{e“^'^'’}). Within this model, the ther¬ 
mal bath consists of N identical subsystems, each of which is 
initialised in the thermal equilibrium state at inverse tempera¬ 
ture /3 for the systems hamiltonian Hg. 

To characterize the bath-system interaction, we assume a 
collision model, where the exchange of information between 
the system and the bath is described by discrete steps. In each 
step, the system S interacts only with a single subsystem of 
the bath. We also assume that the system does not interact 
more than once with the same copy throughout the whole pro¬ 
cess. The one step interaction is described by the map 
acting on the system. The iteration of <i );3 eventually takes the 
system in a state that is arbitrarily close to the thermal equi¬ 
librium state after n steps ll20ll . 

Specihcally, we assume now that is dehned in terms of 
a global partial swapping unitary given by P = cos[0]I + 

i sin[0]r, with 9 s]0,7r/2] and where I = Xg ® is the 
ik) 

identity matrix in Tig ® Til and T is the total swapping uni- 
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tary. Without risk of ambiguity, we can simplify the notation 
by suppressing the index (k). We can now define a thermalisa- 
tion process as a completely positive trace preserving (CPTP) 
map given by <^p{ps) = Tr^{P(ps ® Pi3)P''}- It can be 
easily seen that the local state of the system after the partial 
trace is given by Hill 

^p{Ps) = C^Ps + S^pp+ics[pb,Ps] ( 1 ) 

where c = cos[0] and s = sin[6>]. Note that, in our definition, 
the map represents a legitimate thermalisation model for 
any value of 9 g] 0, 7r/2]. For example, a map correspond¬ 
ing to 9 — t: 12 (i.e. c = 0 and s = 1) would thermalise 
the system in a single step. Clearly, a more realistic model of 
thermalisation requires a vey small value for s. However, our 
results are valid for any value of c and s (we only exclude the 
cases of 0 = 0, i.e. c = 1 and s = 0, since would trivially 
leave any state of the system unchanged). 

It is worth to note that a thermalising machine described by 
the map satisfies two basic requirements of a thermalisa¬ 
tion process; 

1. It is a thermal operation ifl^ . i.e. a CPTP map $(ps) = 

'Yt:b{UPs^ppU^} where (i) the state pp is a Gibbs state, 
i.e. pp = ^ ^ and (ii) [/ is a global unitary such 

that [U,H,'+Hp] =0. 

2. It obeys to zero-law of thermodynamics, in the sense 
that the iteration of the map will eventually bring the 
system to its thermal equilibrium state. 

Zero-law can be formally written as 

Vps,Ve,3no : {n > no} ^ {Di^(^'p(ps)\pp) < e}. (2) 

where represents n subsequent iterations of the map 
It is important to point out that equation dD is an essential 
requirement in our definition of a thermalisation process. 

-Counter-acting thermalisation: The vague concept of act¬ 
ing against thermalisation can be interpreted in different ways. 
Here we consider a precise type of counter-thermalisation pro¬ 
cesses, the restoring processes and the stabilising processes. 
One simple way to intend the concept of contrasting thermali¬ 
sation consists of invoking a restoring process. Such operation 
may be described by a restoring map, i.e. a CPTP-map <I>r that 
allows us to restore the system in a given target state ps at any 
point during the thermalisation process. More precisely, given 
a thermalisation process ^p whose equilibrium state is pp, we 
can define the thermalisation path from ps to pp induced by 
<I >/3 as the set of all the states such that = <I>^(ps) for 
some n. Thus, a restoring map <I>r. is defined as a map that, 
after m iteration, transform any state p} belonging to the ther¬ 
malisation path into the target state ps. 

The conditions defining restoring operations may be con¬ 
sidered quite demanding, in the sense that the map is required 
to restore the state ps when acting on any state in the thermal¬ 
isation path (which include states arbitrarily close to pp). In¬ 
deed, one may consider a wider class of maps, which are only 
able to restore the system after it has gone through a fixed 
number of thermalisation steps. In particular, we consider 


the maps that restore the system only if it has gone through 
a single thermalisation step. Thus we introduce a second class 
of counter-thermalising maps, the class of stabilising maps. 
We call a stabilising map a CPTP-map which ’’continuously” 
keeps a system undergoing a thermalisation process in a given 
target state ps. More precisely, a stabilising map with respect 
to Ps and ^p is a map <I>s : B{'Hs) B{TLs) such that 

$«($p(p,)) =p,. (3) 

Clearly, the stabilising map will depend on the thermalisa¬ 
tion process ^p and on the target state p^. We can very well 
say that, for any thermalisation model and any initial state 
Ps, it exists a class of stabilising maps $s, and we will denote 
this class of maps with ps). Finally, we define the class 

of all stabilising maps with respect to ps, denoted S{ps), as 
the union of all sets 5($p,ps) for all possible thermalising 
machines, i.e. for all possible maps based on the partial swap¬ 
ping operation. 

-Asymmetry and coherence measures: At this stage, a pre¬ 
cise definition of coherence is needed. We define the amount 
of coherence of a density matrix p, with respect to an Hamil¬ 
tonian H, by using the conc ept of sy mmetric operations under 
time-translation. Following lUTHl^ . we introduce the class of 
symmetric operations with respect to an Hamiltonian H. This 
class includes any CPTP-map that commute with the time evo¬ 
lution operator generated by H. That is, a symmetric map <I> 
is such that The set of 

symmetric states with respect to H are then defined as the 
set including any state p that is invariant under time transla¬ 
tions, i.e. e“*^*pe*^* = P- The symmetric states are thus the 
states which are diagonal in the eigenbasis of H (for simplic¬ 
ity, here we assume that the Hamiltonian is non-degenerate. 
All the results can be generalised to the case of a degener¬ 
ate spectrum with no conceptual difficulties). Naturally, we 
will call the states which have non-zero off-diagonal elements 
in the eigenbasis of H asymmetric states. Asymmetric states 
are then states with coherences, and the words asymmetry and 
coherence will be used with the same meaning hereafter. 

The amount of coherence in a state can be quantified by 
any real function C(p) which do not increase under sym¬ 
metric operations, i.e. any function such that VT* symmetric, 
(^(^(p)) < C(p). In particular, here we will use a specific 
coherence measure which is based on the Zi-distance between 
Hermitian matrices H . The Zi-distance is defined as 

Dh (pIp') = X] 1^* J “ P'id I ■ 

i,3 

The amount of coherence Ch (p) in a state p with respect to 
an Hamiltonian H is then given by the Zi-distance between 
p and its projection on the subspace of symmetric states with 
respect to H, i.e. 

Ch{p) = DiMPsip)) (5) 

where Ps{p) = P*.il‘P*)(<Pi|, with {|p>i)} the eigenstates 

of H. Ps{p) represents the diagonal part of p in the basis 
{|<Pi)} (the block-diagonal part in case of degeneracy). Thus, 
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the amount of coherence in p with respect to H is given by 
Ch{p) = \pij\- Notice that, according to this defini¬ 

tion, a state is diagonal (block-diagonal for degenerate Hamil¬ 
tonians) in the eigenbasis of H iff it has zero coherence. 

-Generalised thermal operations: We now introduce the 
class of operations that we are allowed to perform. We restrict 
the class of allowed operations to the class Q of Generalised 
Thermal Operations (GTO) defined as $(ps) = T^rlUps 0 
PrU^} with [U,Hs -f Hr] = 0 and [pr,Hr] = 0, where U 
is a unitary operation and pr is the state of an ancillary sys¬ 
tem R whose Hamiltonian is Hr- This definition represents 
a generalisation of the definition for thermal operations. In¬ 
deed, while we retain the energy conservation condition on 
the unitary U, we now allow the use of ancilla states which 
are not Gibbs states, as long as they are stationary states with 
respect to the local Hamiltonian Hr - Notice that we do not put 
any other restriction on the system R, which represents the 
energy-resource system. This means that we are also allowed 
to use resource systems with arbitrarily big Hamiltonians. 

The motivations behind the choice of considering only 
GTO will become clear shortly. At this stage we only notice 
that the two condition [U, Hg + Hr] = 0 and [pr, Hr] = 0 are, 
if taken individually, neutral, meaning that none of the two 
conditions generates any restriction on the class of allowed 
maps if considered alone. This can be easily proven by tak¬ 
ing U to be the swapping unitary T and then imposing only 
one condition at the time. It is easy to see that in both case a 
map that transforms an arbitrary state p into another arbitrary 
state p' can be easily constructed. This is not true for thermal 
operations. It is also worth to note that GTO are, as well as 
thermal operations, symmetric operations with respect to Hg 
ifr^ . This means that, any proper coherence measure such as 
Ch^ (p) must not increase under GTO. This property plays a 
crucial role in the following derivation. 

-State’s stabilisation and restoring: We are now ready to 
state the main results of this paper. Our first goal is to es¬ 
tablish whether it is possible or not to hnd an operation that 
contrasts thermalisation, meaning either a restoring or a sta¬ 
bilising operation, which is also a GTO. We find that this is 
possible if and only if the target state ps (i.e. the state we 
want to stabilise or restore) has no coherences. Having intro¬ 
duced the class S{pg) of stabilising maps with respect to pg 
and the class Q of generalised thermal operation, we can state 
the following proposition: 

Proposition. A stabilising operation with respect to pg which 
is also a generalised thermal operation exists if and only if 
ChAPs) = 0. Formally 

ChAps) = 0 ^ s{pg) r\gf% ( 6 ) 

Proof Here, we give a sketch of the proof. The full details of 
the proof are founded in appendixiAl 

In order to prove the left-to-right implication (i.e. zero co¬ 
herence is a sufficient condition for the existence of the map) 
we explicitly construct a stabilising map which is also a GTO. 
Such map is constructed by taking Hr = Hg, pr = Ps and the 
swapping operation as the global unitary U. 

In order to prove the right-to-left implication (i.e. that zero 
coherence is a necessary condition for the existence of the 


map), first we prove that, for any ps with non-zero coher¬ 
ence, strictly decreases the coherence of pg. Thus, a sta¬ 
bilising map must necessarily increase the coherence of 
^p{ps). On the other hand, it is proven that GTO cannot in¬ 
crease the amount of coherence in any state lfl^ . Hence, it 
follows that a stabilising map that is also a GTO cannot exist 
when ChAPs) ¥= 0- □ 

Clearly, this result is extended to restoring maps too, since 
these maps are also stabilising maps. Indeed, in the case of 
restoring maps the result is even more general: it can be easily 
proven that a restoring operation which is also a GTO can¬ 
not exist not only for thermalising machine, but for any model 
that obey zero-law of thermodynamics as stated in equation 
dU if the target state pg has non-zero coherences. This can 
be easily understood by noticing that such operations must be 
able to increase coherences in order to transform states which 
are arbitrarily close to pp into pg . Since GTO cannot increase 
the amount of coherences, it follows that a restoring map <I>r 
cannot be a GTO. 

We now pass to analyse the consequences of proposition 
dUl. It is clear that, if we are willing to maintain a state with 
coherences out of thermal equilibrium, we need to abandon 
one of the two conditions defining GTO. Since the resource 
system R has the role of an energy reservoir, it is natural 
to keep the condition [pr^ Hr] — 0, which ensures a precise 
knowledge of the amount of energy available in the reservoir 
at the initial time. This means that we need to consider uni¬ 
tary operations that do not commute with Hg -f Hr ■ We thus 
assume that in addition to the energy available in the resource 
system, the operator is able to change the total energy by tun¬ 
ing the coupling between the resource system R and the sys¬ 
tem S. This energy can be interpreted as work performed by 
the operator on the system {S -f i?} J^^-lUt]. 

Let us then consider maps 'hs(Ps) = Tir{Upg ^ PrUA for 
which [pr, Hr] — 0 but any unitary U is allowed. To simplify 
the notation, here we have set p'^ = <I>; 3 (p;,). The work W 
done by U on the Sh-R composed system can be written as 

W = Tr{(Hg + Hr)Up', ® prUA - Tr{(iT. + Hr)p'g ® Pr} 

(7) 

The connection between W and the commutativity between 
Hg -f Hr and U is then clear: proposition (|6]l ensures that, 
when the state pg has zero coherence, it is always possible 
to find a U that commute with Hg -\- Hr, in which case we 
have kk = 0, as can be easily seen from equation (|7]). On the 
other hand, proposition (|6ll also ensures that we must choose 
an unitary such that [Hg Hr, U] A 0 in order to stabilise 
a state with coherences (assuming that we keep the condition 
[pr. Hr] = 0), in which case we have in general that W f Q. 
The quantity W can be thus interpreted as the work needed in 
order to preserve coherences from thermalisation. 

To further develop this concept, we now present an example 
of a map which stabilise states with coherences and we 
explicitly calculate the value of W. The stabilising map 
is constructed by considering a resource system composed by 
many copies of the system itself, all prepared in the target 
state Pg and a global unitary given by the swapping operation 
T. The map is then given by $s(Ps) = Trr{Tp[. 0 Ps^A 
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with T the total swapping operator, whose action is given by 
Tp'g® PsT^ — Ps® p's - In order to ensure that [ps,H,r\ — 0, 
we pick an Hamiltonian H^- for the resource system given by 

Hr = -(l//3)(log[^Hj + log[ps])- (8) 

with ZHr = the partition function defined by 

Hr and the temperature /3. Clearly, Hr is constructed such 
that ps = is the corresponding Gibbs state at 

temperature /3. We can now calculate how much extra work 
W is required to perform the stabilisation. Substituting equa¬ 
tions o and ([8]) in equation Q and writing Hg as Hg = 
— (l/^)(log[Z//J -I- log[/ 0 / 3 ]), it can be shown with a few 
straightforward passages (see Appendix iBli that the work W 
can be written as 

W = — L>symm(Ps 1/0/3 ) (9) 

where f9symm('|') is the symmetrized relative entropy de¬ 
fined as £>symm(p|p') = D{p\p') + D{p'\p), with D{p\p') = 
Trjplogp — plogp'} the usual relative entropy. Notice that 
the quantity in Equation @ is always positive, which means 
that external work is always required by the system in order 
to maintain coherences, as expected. 

-Conclusions: We have shown that, given a thermalisation 
model based on partial swapping operations and described by 


the map a stabilising map which counteract and 
which is, at the same time, a generalised thermal operation 
exists if and only if the target state ps is block diagonal in the 
system’s Hamiltonian eigenbasis. This implies that, in order 
to keep a state with coherences out of thermal equilibrium, 
one of the two conditions defining the class of GTO must 
be violated. In this regard, we have analysed maps such that 
\U,Hg -\- Hr] ^ 0, finding that an extra amount of work W 
proportional to the symmetrized relative entropy between ps 
and pp is required to perform such operations. The work W 
can be then naturally interpreted as the amount of work that 
is necessary to perform on the system in order to maintain the 
coherences in pg. This result, albeit derived in the context of 
collision models, can be extended to the case of a continuous 
process in which a system is allowed to interact with a bath 
and a resource system at the same time. Indeed, a continuos 
process can be simulated by a collision model in the limits of 
infinitesimal steps lEhU^ . 
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Appendix A: Full proof of Proposition 

Proposition. A stabilising operation with respect to ps which is also a generalised thermal operation exists ijf ChSPs) = 0- 
Formally 


{Ch. {Ps) = 0} ^ {Sips) n a ^ 0} (Al) 

Proof. We first prove the left-to-right implication (Sufficient condition) and then the right-to-left implication (necessary condi¬ 
tion). 

-Sufficient condition: First we prove that Ch^ (Ps) = 0 is a sufficient condition for the existence of the map, i.e. 

Gan5(p«)}. 

To do that we explicitly construct a map $s(Ps) which take any state p'^ (including ^p{ps)) into ps and that is a generalised 
thermal operation. We use an ancillary resource system R whose initial state is pr and whose Hamiltonian is Hr. The map 
$s(Ps) '^hen defined as 

‘f.(p;)=Tr,{rp'®p,Tt}, 

where the unitary T is the total swapping operation whose action is given by 

Tp'g (g) PrT^ = pr® p's- 


By choosing pr = Ps and Hr = iFs, is a then a stabilising operation with respect to ps and it can be easily checked that it 
is also a generalised thermal operation (the swapping operation is energy preserving since Hr = Hg and [ps, Hg] = 0 since the 
state Ps is bock diagonal in the eigenbasis of Hg). 

-Necessary condition: we now prove that {ps) =0 is also a necessary condition for the existence of a stabilising map, i.e. 

Step 1: As a first step we prove that for any map written in the form $; 3 (ps) = Tr; 3 {P(ps ® pis)P'^}, with P = 
cos[0]I -f i sin[0]T and 9 g] 0, 7r/2], and for any state ps with non-zero coherence, we have that 

CnA^piPs)) < ChAPs)- (A2) 

To do so, we explicitly write the state ^ffipg), which according to equation ([T]i is given by 

^p{ps) = C^Ps -f S^pp -f ics[pp, Ps]. 

We now pick the basis {|t)} of eigenvectors of Hg in which pp is diagonal (although here we consider non-degenerate Hamil¬ 
tonians for the sake of simplicity, choosing such basis is necessary to extend the prove to the case in which Hg has degenerate 
eigenvalues). Following the definition of coherence given in Equation (|5]l (i.e. the sum of the modulus of the off-diagonal 
elements in p) we can easily derive that 

CnA^piPs)) = I (A3) 

where p^'^d) = (i|p|j) are the matrix elements of the two density matrices in the basis {|i)}. Notice that, since represents 

probabilities, we have that — p^p'^'^A ^ 1- Taking into account that we consider maps with 9 G]0,7r/2], which means 

c G [0,1[ and s g] 0, 1], we can write 

CnA^piPs)) = Y1 lAp^g'"'’ +ics{p^p’''^ - Pp'^^)pA^'’\ 

= Y1 \ cPs’^^\\c + isiPp’'^ -Pp’^^)\ 

= ^^\Ps’^^\\/^'^ + (Pp - Pf^^f < ^C|p^^)|v^c2-f s2 

=cY,\pA^'>\ = c-ChAPs) 

iAj 


(A4) 
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Since c < 1 for any non trivial process, we have that 

Ch^{^p{Ps)) < c ■ ChAPs) < ChAPs) 


(A5) 


This prove equation (IA21i . 

Step 2: As a second step, we notice that generalised thermal operation cannot increase the amount of Coherence in any state 
Ps, i.e. for any generalised thermal operation we have 

ChA^s{Ps)) < ChAPs) (A6) 

This property is proven in IT^ . 

Step 3: The last step of the proof is straightforward. Indeed, a stabilising map is such that <i>;,(<i>;3(/9s)) = Ps- Because of (IA21 i. 
must strictly decrease the coherence for any state. Therefore, the stabilising map must necessarily increase the amount of 
coherence in <i>; 3 (ps). Because of (IA61 l. this means that it cannot be a generalised thermal operation. 

□ 


Appendix B; Derivation of Equation @ 


We start by considering the expression in equation (|2l). Keeping in mind that we have 


Pr — Pst 

p's = ^/3(ps) = C^Ps + S^Pfj+icslpfj.pA 

U = T with T representing the swapping unitary 

Hs = -(l//3)(log[ZrrJ + log[p;3]) with Zh, = Tr«{e“^^'’}, 
Hr = -(l//3)(log[^rrJ + log[ps]) with Zur = Trs{e“^^'-}, 

equation (|7]) can be transformed according to 


W =Tr{(H, + Hr)Up', (8) PrU A- Tr{(Tf« + Hr)p, ® p.} = 

Tr{(iT, + Hr)TA^ ® PsTA - Tr{(Tf« + Hr)p', ® pj 
=Tr{(iTs + Hr)ps ® pA - Tr{(7Ts + Hr)p'g O pj 
=TVs{HAPs - p's)} + TVriHrip's - Ps)} 

=TVs{Hsips - Aps - S^pp - ics[pp, Ps])} + TXr{Hr{Aps + pp + ZCs[pp, Ps] - Ps)} 


= 5 


Trs{Tfs(Ps - Pp)} - TXriHrips - pp)} 


— ICS 


Trs{[Tfs,Pp]Ps} + TTr{pp[ps,Hr]} 


(Bl) 


(B2) 


Noticing that [Hs,pp] = [Hr,ps] = 0, the second term in the last line of equation (IB2I) vanishes. After substituting the 
expressions for Hg and Hr in the first term, we have 



logZ//^Trs{(ps 


Pp)} + log ZnAttrAPs - Pp)} - Trs{(ps - pp) logpp} + Tr^Kps - pp) logps} 


(B3) 


Since ps and pp have are both trace-1 matrices, the first two terms in the expression above vanish. Re-arranging the two 
remaining terms and considering that the system r is a copy of the system s, so that tracing over s and tracing over r are 
equivalent, we obtain 


r 

W =—[^-Trs{(ps - pp) logpp}-fTr^{(ps -pp)logps} = 
= ^ [Tr{ps log Ps - Ps log pp} + Tr{pp log pp - pp log ps} 
= ^[D{Ps\pp) + D{pp\ps)]. 


(B4) 


This is the expression given in equation (|9l) for W. 











